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Abstract 

We propose a new treatment of the Hubbard model that is based both on the coherent-potential 
approximation (CPA) and the virtual-crystal approximation (VCA). It is well known that the equi- 
librium found using the one-particle CPA Green's functions does not predict an ordered magnetic 
ground state, while Stoner's mean-field treatment, which is equivalent to the VCA on the Hubbard 
model, does so for a wide range of parameters. A hybrid treatment, the r-CPA, is developed, in 
which a particle is assumed to be scattered from an array of static opposite spins for a time r 
related to the inverse of the band width. The propagation is treated in the CPA over this period; 
thereafter the particle sees the time-averaged effect of the scatterers and hence can be treated in 
the VCA. This model, with suitable approximations, does predict magnetism for a modified Stoner 
criterion. 

PACS numbers: 71.10.Fd,75.10.-b,75.40.Cx 
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I. INTRODUCTION 



While much progress has been made in recent years in the theoretical treatment of mag- 
netism in the transition metals and similar materials, a complete treatment of the many 
aspects of these familiar but complicated systems is still lacking. The electrons retain many 
of their energy band properties as is reflected by their metallic behaviour, but the mag- 
netic properties reveal that strong correlations exist between the electrons on the individual 
atoms, and exchange occurs between them. 

The simplest treatment of metallic ferromagnets is due to Stoner— who assumed that 
the electron interaction gave rise to a molecular field which predicted ferromagnetism in a 
simple band model if the ratio of the interaction to the band width (or more precisely the 
density of states of the Fermi surface) was adequately large. But beyond the mean field the 
model gives no correlation between the magnetic electrons. More complex band structures 
also predict antiferromagnetism in some cases- 

In recent years much more sophisticated treatments of the band structure using estab- 
lished methods such as those based on the spin dependent local density approximatio n 3 ' 4 ^ 6 
have been able to make much better predictions of the Fermi surface and related magnetic 
properties, but again a full treatment of the atomic electron correlations is lacking. Hybrid 
treatments such as that of Gutzwilleri^ which start from a mixture of correlated and un- 
corrected states can give an improved treatment of the magnetic properties at the expense 
of some of the metallic ones. 

An alternative approach to this problem has been through the study of simple band mod- 
els with electron interactions of which that due to Hubbard^ is the simplest which appears 
to contain the essence of the problem. He restricts the electron-electron interaction to occur 
on single atomic sites which in the case of a single band only occurs between electrons of 
opposite spin because of the exclusion principle. Hubbard, in a series of paper o 10 ' 11 , pio- 
neered different approximate treatments to this model Hamiltonian. The simplest replaces 
the electron-electron interaction by its average so that an electron of one spin sees a field 
proportional to the density of the other spin and the system has simple energy bands sepa- 
rated by an energy proportional to the magnetisation and hence is completely equivalent to 
Stoner's mean field approximation. In a different vocabulary it can be regarded as a virtual 
crystal approximation (VCA) where the electrons move in energy bands derived from an 
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average potential. 

In order to give a better treatment of the electron-electron scattering Hubbard^ pioneered 
another approximation where it was assumed that electrons of one spin were frozen on their 
sites and the electrons of the other spin moved through the random potential so created. He 
proceeded to solve this problem by determining an effective potential (which was complex 
and energy dependent) which took account of the average scattering from this array in an 
analogy with an alloy. This method has been independently used in a wide variety of alloy 
problems and is usually referred to as the coherent potential approximation (CPA) 12 ' 13 ' 14 ' 15 . 
It is clear that this approximation over-estimates the electron-electron scattering and in 
particular causes the energy bands of each spin character to split if the interaction is large 
enough. As a result the tendency to ferromagnetism is under-estimated and it has been 
demonstrate d 16 ! 17 that no magnetic order is predicted except in the special case of a half 
filled band. 

Several attempts have been made to improve on these approximate treatments, for ex- 
ample using a dynamical CPA^S^ based on the integration of this approximation into func- 
tional integral techniques. There has been a number of other attempts to modify the alloy 
analog y 20 ' 21122123 . In this work we propose a simple physical model which is intermediate 
between the straightforward VCA and CPA treatments. In it we assume that the random 
scattering centres experienced by an electron of one type in the alloy analogy exist only for 
a certain time r which is proportional to a typical hopping time between sites and hence 
inversely to the band width. At later times the electrons see a completely average potential 
reflecting the movement of the scattering centres. Thus the propagators are calculated in 
the CPA approximation for the short time interval r and have to be matched to the VCA 
propagators which hold for longer times. Although this is a straightforward procedure in 
principle certain further approximations are necessary to make it tractable. The properties 
of the system are evaluated in this so-called r-CPA. The single particle properties such as 
the density of states show, as expected, a mixture of characteristics associated with the 
VCA and CPA. Ferromagnetism is predicted to occur for a range of parameters which are 
somewhat more restrictive than those in the simple Stoner model. 

This paper is organised as follows: Sec|n] introduces the model used. The VCA and CPA 
methods are outlined in Sec lIIII The r-CPA, mixing the VCA and CPA, is constructed in 
Sec JIVI Numerical results for the density of states and the magnetisation using a simple 
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cubic band are presented and discussed in SecfVl 



II. MODEL HAMILTON! AN 

The Hubbard model is currently widely considered to be the most concise model that 
captures the essence of the interactions in an electronic system. It allows the hopping of 
the electrons from site to site to compete with the Coulomb repulsion while including the 
Pauli principle. We will consider the simplest case of the Hubbard model, a single-band 
with nearest-neighbour hopping and on-site repulsion : 

H = Yl ^ c ta c ^ + U Yl c i^T c S c a (!) 

i,j,cr i 




cj a and Cj CT are the usual creation and destruction operators for electrons in the Wannier 
states. The free part of the Hamiltonian H is evaluated in the Bloch states H = J2k e k c tcr c ka 
using a tight-binding approximation for the band structure e^. We will concentrate in this 
paper on the simple cubic band of half band width w — 1. As for the interaction, it is 
repulsive and only felt when one electron of spin a meets another of spin — a at the same 
site. This model simulates relatively well a real material possessing narrow energy bands, 
like the transition metals and their alloys. We will work on a three dimensional lattice of 
N = L 3 sites, at T = 0. We use the framework of the canonical ensemble and fix the number 
N e of electrons per atom that distribute themselves at random over the lattice sites. The 
Fermi statistics takes care that no site is occupied by two electrons of the same spin while 
the interaction favours singly occupied sites. The propagation of an electron of spin a in 
the lattice can be described by the retarded Green's function 

G-(M') = - l 2n6(t-t>)({cUt),c+ a (f)}+) 

=:«c*(f);c£(0» (2) 

where H has been set to 1. The following Fourier transform (where no particular origin is 
specified) 

1 r°° 

G a (E) = — d(t- G a {t - t')e iE ^ (3) 

27r J-oo 

will often be used, and t' will be set to 0. 

We are interested in the properties of this system averaged over all possible configurations, 
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for which the Green's function G^(E) can be defined by 

k 

The density of states p a is obtained from the single-particle Green's function with 

p°(E) = —MF a (E) (5) 

7T 

where F a (E) is the trace of the Green's function. It can be obtained as 

F ^ E ) = ^Y, G ^ E ) (6) 

k 

The relative proportions of the number of up-spins and down-spins must be determined 
self-consistently together with the chemical potential fi. Defining and respectively 
as the average number of spins up and down per atom, the equilibrium conditions read 

N a = ^ dEp a (E) (7) 

J — oo 

with a E {T, 1} 



III. VCA AND CPA 



The virtual-crystal approximation and the coherent-potential approximation have been 
both applied with various success to both the Hubbard model and the somewhat related 
problem of the random binary alloy A\- C B C . It is worth recalling briefly the development of 
the VCA and the CPA in the binary alloy before turning to the Hubbard model. 



A. VCA and CPA in the binary alloy 

The analogy between the Hubbard model and a binary alloy was noted by Hubbard in 
the last paper of the series^. We consider a particle propagating on a lattice where each site 
is occupied at random either by a host atom A or an impurity atom B. The concentration 
of impurities is c. The model is represented by H a u oy = H Q + V, where V is a quadratic 
operator of the form V = Ylm \ m ) e «* ( m \- The atomic potentials can take two values, either 
€a or €b, depending on whether the particle is on an A site or a B site. Setting to and 
e B to U, we find that this is very similar to the Hubbard model (hence the alloy analogy), 
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for which the spin a propagating in the lattice interacts with energy U only on sites where 
a —a spin is present. The analogy would be exact if the spins of opposite direction could be 
considered as frozen in the Hubbard model. 

The VCA and CPA are well established method o 15 ' 24 and will only be outlined here. 

If P describes the propagator in the perfect crystal defined by H , the perturbation expansion 

in the whole Hamiltonian H a u oy reads 

G = P + PVP + PVPVP + --- = P + PVG (8) 

The problem of evaluating the configurational average (G) = P + P(VG) is solved in the 
VCA by the random phase approximation decoupling (VG) ~ (V)(G). With (e m ) = cU, 
the propagator for the binary alloy in the VCA is 

Gl CA {E) = 1 (9) 
hi — ek — cU 

While the VCA is only a valid approximation for small perturbations V, the CPA is a par- 
ticularly successful method that interpolates between the limits of strong and weak disorder 
and interaction. The CPA assumes the existence of an effective medium described by the 
Green's function G e and self-energy £ that obey the Dyson equation 

G e = P + PSG e (10) 

The true Green's function G can be developed in the perturbation expansion (jSJ), which, 
together with f(10j) . lead to a new expansion of G in terms of the medium 

G = G e + G e {V - E) G (11) 

A particle moving in the effective medium sees the impurities as embedded in a uniform 



potential. As shown in Refs ll5ll24 it is also possible to expand the true G in terms of G e 



and a scattering T matrix so that 

G = G e + G e TG e (12) 

The T matrix can be expressed as a function of the individual atomic scattering matrices 
T 



- Ill 



T — T m + T n G e T m + T n G e T m G e Ti + ■ ■ ■ (13) 
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Projected on the the Wannier states, the scattering on the m atom can be found to b 



,15.24 



tm ~ 1 - (7 m - S)F (M) 

F is the diagonal element of G e , F = (m \ G e \m) (see EJ). A central step of the CPA is to 
demand that the medium propagator G e corresponds to the true propagator averaged over 
the disorder: G CPA := (G) = G e . Consequently, from (fT2j) . the total scattering of impurities 
is set to zero: (T) = 0. (T) can then be decoupled by a random phase approximation on 
the single T m matrices. We note that this is a much better approximation than the VCA, 
since the repeated scattering by the same site is forbidden, as (JT3j) shows. In this form, the 
CPA neglects the scattering from clusters of atoms, though some extensions of the theory 
do include this effect. 
The CPA condition becomes therefore 

(tm) = (15) 

This is the central equation of the CPA. The self-energy is calculated such that the scattering 
generated by any single impurity is zero on average. The CPA approximation for the one 
particle Green's function is then given by 

with the self-energy £ determined self-consistently with F from (|15jl and p4jl . From these 
two equations and for the binary allow with = 0, e_e = U, the relations between F and E 
can be expressed as 

£= !-(,/%*■ <17) 
where F is the diagonal element of G^ PA . While the VCA self-energy is a real number that 

merely shifts the energy levels by ell, the CPA S is an energy-dependant, complex number 

which accounts for the damping of the quasiparticle states by the impurities. 

B. VCA and CPA in the Hubbard model 

The VCA for the Green's function (j2J) in the Hubbard model $IJ can be obtained in two 
different ways. The first approach consists in applying the Hartree-Fock approximation to 
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the Hubbard model and effectively reducing it to a Stoner model 

Hstoner = H + U ( N ^H + ( 18 ) 
i 

The usual number operator notation h ia = cf a Ci a has been used. The Stoner model (fT8j) 
leads directly to the set of Green's function analogous to © 

GZ VCA (E) = (19) 

The up-spins see the (frozen) down-spins as impurities of concentration with whom they 
interact with energy U, and vice-versa for the down-spins. 

The second way of obtaining the virtual-crystal approximation is to decouple the equation 
of motion for the Green's function (j2J) in the Hubbard model 

d 
dt 



i-fA(Cia(t),Cp) = 27T<y(*)<5y + ^ t U ( (c ltT (t) , C+ ) ) 



i 

+ U{(h i ^ a {t)c ia {t),cj a )) (20) 

If h,i_ a (t) is taken out of the last Green's function in (J2Uj) and averaged separately, 
(fii-oif)) ~^ N~ a , the equation of motion is decoupled and the solution is given again 
by (|19j). If, on the other hand, ni_ a (t) is considered at t — and takes the value or 1 at 
random, applying thus the alloy analogy, the CPA result (|16I17|) can be retrieved. 

1 



Gl CPA {E) 



S ^ = 1 - [U - Z°(E)}E°(E) ^ 

This is however not a straightforward procedure. Hubbard developed an elaborate scheme 
based on the alloy analogy in the last paper of the series^, calling it the "scattering correc- 
tion". It was noted later by Velicky et al that the scattering correction Green's function is 
the same as the CPA result for the alloy model. 



IV. THE r-CPA 



It has been known for some time now that the CPA does not allow for ferromagnetic 
solution a 16 ! 17 . On the other hand, the VCA treatment of the Hubbard model leads to 
an artificially strong ferromagnetic phase. We propose to combine these two different 
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approximations in order to keep the advantages of the CPA while allowing for spontaneous 
magnetic ordering to occur. As was seen in the previous section, the VCA and CPA take 
two very different views when considering the motion of a particle in the lattice. The CPA 
sees the scatterers as fixed at the positions given at t — 0, while the VCA considers the 
scattering events only as averaged over a time t — > oo. In other words, the CPA is valid for 
the Hubbard model only at t = and the VCA at t — > oo. A better physical picture should 
emerge if the CPA treatment is applied for a short time only, when it is still reasonable to 
consider the scatterers as fixed. This time is typically the time needed for the particle to 
travel between sites. Thereafter we can assume that the VCA average provides a reasonably 
good treatment. 

Let G\ CA {t) and G k PA (t) be respectively the time Fourier transform of the VCA and 
CPA Green's functions for t > 0. The central assumption is that up to a time r > and 
t « 1/w the system is described by the CPA Green's function, whereas beyond r it is 
described by the VCA Green's function. A new function gl CPA (t) is defined as 

9l CPA (t) = -i2n 9{r - t)G c k PA (t) - a *2tt 6{t - r)G v k CA {t) (22) 

At t = 0, g k CPA is equal to the CPA Green's function and is thus properly normalised. A 
factor a, the "matching factor", has been introduced to ensure the continuity of the wave 
function at time r. 

The retarded Green's functions for the VCA and CPA are given respectively by G% (E) 
(J!2J) and G k PA (E) (JTHJ). The Fourier transform reads: 

/oo 
dEe- iEt G k {E) (23) 
-oo 

It is straight-forward to obtain G\ CA {t) by contour integration in the lower half of the 
complex plane. There is one pole at z = e k + cU — iS, so that 

G v k CA (t) = -t2n 6(t)e- i{ek+cU)t (24) 

It is more difficult to do a similar calculation for the CPA. The problem lies in the nature 
of the CPA poles. It has been noted by Velicky et that the two complex CPA poles 
cannot in all generality be interpreted as quasiparticle energy. A calculation for the CPA 
similar to the one made for the VCA does not lead to a correct result. In particular, the 



correct inverse Fourier transform is not retrieved. In order to keep the denominator of the 
energy Green's function as E — e& — we assume instead, by similarity with the VCA: 

Gf PA (t) = -i2it 6(t)e~ i{tk+s)t (25) 

At this stage the fact that S is an energy-dependant quantity is ignored, and S is treated 
as a fixed quantity rather than a function. By doing so, Fourier-transforming (J25|) gives 
the expected result JTEJ). The substitution £ — > S(-E) must be made thereafter. The whole 
gl CPA (t) is still normalised at t = 0, but because of the dependence of £ on the energy, 
this does not guarantee the normalisation of gJ, CPA {E). The normalisation will have to be 
corrected numerically. 

The integrals for the Fourier transform on the partial intervals give, without taking into 
account the matching factor a, 

QT\E) = -L dt e^G\ c \t) = | (26) 

2tt J E - e k - cU 

1 r°° 1 _ p i{E-e k -T,)r 

S C k P \E) = i- jf dte^ PA (t) = \_ ek _^ (27) 

The Q\ CA {E) and Q^ PA {E) are the standard VCA and CPA Green's functions, weighted 
by a complex r- dependant quantity. 

The factor a that forces the CPA and VCA Green's functions to match at t = t is found to 
be 

a = e-^- cU)T (28) 

This value for a can be reintroduced into (|22j). The result for gl CPA (E) is then obtained 
from the Fourier transform of 



gf PA (E) = ±- J Q dtGf PA (ty Et 

1 _ pi{E-e k -T,)r i{E-e k -T,)r 

= — + (29) 

S is now read as S(£ l ), as calculated by the pure CPA equation (|17j). The expression (|23|) 
now reproduces exactly the CPA behaviour in the limit r — > oo. The VCA approximation 
is retrieved in the limit r — > 0. 



We recall that the normalisation of the Green's function has been lost by choosing to 
ignore the energy dependence of S. It is therefore necessary to calculate numerically for 

10 



each k the factor so that (3^ 1 gJ, CPA (E) is normalised. 

The result ()29)1 is however not giving a physical density of states: p(E) changes sign along 
the E axis. By stopping the Green's function at a time r we have in fact introduced a 
non-physical oscillatory term coming from the factor weighing VCA and CPA 

The oscillatory term destroys the analytical properties of the CPA and VCA. The imaginary 
part of A -1 —ir~ 1 gl CPA (E+iS) does not give the physical density of states, as the branch 
cut no longer runs along the x-axis. It is thus necessary to smooth out the Green's function 
by putting the oscillatory term to unity. We obtain finally 

1 / 1 — e I»E(B)r e InE(£)r \ 

GlCPA(E) = Tk \E-e k -E(E) + E-e k -cU) (31) 

This function now possesses all the properties required. It interpolates between a pure VCA 
Green's function at r — > and the CPA Green's function at r — > oo. We observe that (|3*T|) 
is decomposed into a CPA and VCA contribution weighted by a real r-dependant factor. 
Although r is not strictly determined by the problem, a reasonable value for it is given on 
physical ground by the band width. 



V. NUMERICAL RESULTS AND DISCUSSION 

In this section, the r-CPA Green's function (13 lj) is used to obtain the density of states for 
both species of spins in the Hubbard model. Various values of the impurity concentration 
c (c corresponds to either or N^), of the interaction U and of r are considered. Cal- 
culations are made for a simple cubic band with w = 1. Another one-particle property is 
then evaluated, the magnetisation m versus the interaction strength U, for various range of 
parameters. 



A. Density of states in the r-CPA 

The density of states is obtained from the Green's function (|3*Tjl using (|5I6|) . Fig|T] shows 
the density of states obtained for a strong disorder and interaction strength, U = 6.0 and 
c = 0.4. t is varied from (pure VCA) to r — ► oo (pure CPA). At r = the spectrum 
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FIG. 1: DOS p for the r-CPA using a simple cubic band for U = 6, c = 0.4 and w = 1. 

U=0.8, c=0.3 

0.8 
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FIG. 2: DOS p for the r-CPA using a simple cubic band for U = 0.8, c = 0.3 and w = 1. 

exhibits the single VCA band, centred at = c£7. As r increases, the two CPA sub-bands 
appear, at E = for the main band and E = U for the impurity band. Weight is transfered 
from the VCA band to both the CPA sub-bands, affecting all the three sub-band heights, 
but not their general shape. At r — > oo the VCA band has disappeared and only the two 
CPA sub-bands survive, with their original, pure CPA respective weight. In this graph the 
r-CPA merges with the pure CPA at r w 600. 

More dramatic effects on the band shapes occur at smaller U, when the bands mix. Fig|2l 
for U = 0.8 and c = 0.3, depicts a case for which U and c are sufficiently small so that the 
pure CPA sub-bands are not split. As r increases, the height of the VCA band decreases 
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U=2.0, c=0.2 










x=0.5 




T=1.0 




x=2.0 




1=00 




(a) small values of r 



U=2.0, c=0.2 




(b) large values of r 

FIG. 3: DOS p for the r-CPA using a simple cubic band for U = 2, c = 0.2 and w = 1. 

while that of the main CPA sub-band increases. The general picture produced in that case is 
that for r > 0, the degenerate states near the centre of the VCA band start to scatter, while 
the edges of the band are still untouched. For larger r, the states continue to diffuse to the 
impurity band from the centre to the right edge. At the same time, states of lower energy 
are being occupied, re-centring the created sub-band around 0. Cases at intermediate U are 
not so smooth, as can be seen in Figs |3(a)1 and |3(bj| where U = 2.0 and c = 0.2. In this 
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figure the density of states of the r-CPA becomes equal to that of the CPA for r « 250. The 
mixing of the bands produces irregular features as the diminishing VCA meets the increasing 
main CPA band at larger r, Fig |3(b)| This is due to the fact that a large imaginary part at 
an energy E favours the CPA behaviour compared to that of the VCA. The CPA self-energy 
has an imaginary part only within the CPA bands. When the VCA band is centred within 
the main CPA band but is sufficiently distinct from it, its flat top can be sharpened if close 
to the top of the CPA band, producing the small peaks in Fig |3(b)[ Such drops in the VCA 
part of the density of states also occur for the parameters of Fig|21 However under these 
conditions the VCA band is only slightly shifted from the CPA band, and the sharp peaks 
are compensated by the CPA part. 

B. Magnetisation in the r-CPA 

The r-CPA offers a treatment that leads to density of states that interpolate between the 
pure VCA and the pure CPA. As will be shown in this section, the method also interpolates 
the value of the magnetisation versus the interaction strength U between the VCA ground 
state (ordered ground state for sufficiently large U) and the CPA ground state (no ordered 
ground state at any U). We only consider in this paper the competition between paramag- 
netic and ferromagnetic ground states. 

In practice, the equilibrium conditions for a certain U and a given, fixed number of spins 
N e , are found by iteration. For example, a value for is given as an input. This al- 
lows the determination of (J3*T|) and the density of states from ()5I6|) . N^-, deduced from 
N e = + N^, determines the chemical potential [i from (J7j) with o =|. also gener- 
ates G>, which together with the chemical potential, determine a new value for N\ The 
procedure is repeated until converges. The paramagnetic solution = is always a 
solution. If another solution is found, its stability must be tested against the paramagnetic 
case by measuring the total energy of both solutions. The total energy E of the system is 
found by summing up the separate energies of the two spin a 16 ' 20 



where A a (k,u) is the spectral density defined by A a (k,u) = — n 1 MG^,(uj). Diagrams of 
magnetisation m = — versus interaction strength U can be calculated for any value 




(32) 
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of N e . In Fig0] the magnetisation is depicted for a range of values of U. r is fixed to 1 
and the curves calculated for various filling. The first observation is that at no value of N e 
(N e < 1) is there saturated ferromagnetism. This is due to the residual CPA part around 
E = that is never totally empty. We can also see that at N e pa 0.4 the maximum of the 
magnetisation is not at large U: m increases at around U ~ 8 to reach its peak at U pa 2. 
Keeping the filling constant (iV e = 0.4 and N e = 0.7 respectively), while varying r give the 
plots of FigsEl and El The magnetisation is interpolated between its VCA value at r = 
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Ne=0.7 




FIG. 6: Magnetisation m versus interaction U for the rCPA at various r but a fixed filling of 
N e = 0.7. 

and the CPA value at r = oo. At r = and for the fillings considered, the system is fully 
polarised at large enough U. The magnetisation has already dropped to zero for any U 
at r = 10 for N e = 0.4, and at r = 50 for N e = 0.7. The increase in m before the drop 
in magnetisation is particularly acute for N e = 0.4 and r = 2.0, and is an artifact of the 
model caused by the position of the chemical potential in the VCA and CPA sub-bands. 
Fig0 shows the respective positions of the CPA and the VCA sub-bands and the chemical 
potential for N e = 0.4, compared with N e = 0.7 at r = 2.0. At large U, the density of 
states of the up-spins (the minority spins) exhibits three well separated sub-bands for both 
fillings N e = 0.4 and 0.7. The lower CPA sub-band is completely full. As U decreases, the 
VCA band gets closer to the lower CPA sub-band, in both the up- and down-spin DOS. 
At first, only the down-spin bands are affected, meaning that the chemical potential travels 
down the energy scale while the proportion of up-spins and down-spins remains more or 
less constant. This is the situation illustrated in the two uppermost plots of Figd When 
U is further decreased, the chemical potential will eventually reach the edge of the CPA 
lower sub-band. At a filling of N e = 0.4, this is happening at a fairly large U, U ~ 8.0, so 
that the VCA band is still well separated (plot in the middle of the left-hand side, FigEJ). 
N] will thus decrease until the VCA band merges with the CPA band (bottom plot on the 
left, FigEJ). For a higher filling, the two sub-bands merge before or at the same time as the 
chemical potential gets into the CPA band. The density of states at the chemical potential 
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FIG. 7: Density of states for both the up and down particles in the rCPA with r = 2.0. On the 
left the concentration is kept at N e = 0.4 for U = 10.0, 8.0 and 3.0. On the right N e = 0.7 and 
U = 8.0, 2.8 and 2.3. Occupied states are marked by the hatched area. 

thus rapidly increases for the up-spins, and the magnetisation decreases accordingly as the 
up-spin band matches the down-spin band. 

The increase in m at intermediate U, though also occurring at other small fillings, becomes 
dramatic around N e m 0.4 and r m 2.0. It is at this order of magnitude of r that the system, 
at iV e ~ 0.4 , exhibits a balanced CPA-VCA behaviour. At much smaller r, the relative 
weight of the CPA band becomes very small compared to the VCA part. The chemical 
potential falls within the CPA sub-band of the up-spin over a large range of U before the 
up-spin VCA band reaches it. At large r the system becomes more CPA in character with 
a small VCA part responsible for the low magnetisation. 

VI. CONCLUSION 

A new treatment of the Hubbard model is proposed that combines the physical pictures 
and respective advantages of the CPA and the VCA. It is based on the fact that the CPA 
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on the Hubbard model considers the scatterers as fixed at their positions at t = 0, whereas 
the VCA performs a time-average of the system. The one-particle Green's function in the 
so-called r-CPA is built from the one-particle time Green's functions of the CPA and VCA. 
The CPA Green's function is applied up to a time r, thereafter the VCA Green's function, 
multiplied by a factor for continuity, is used. This time r is taken to be of the order of 
magnitude of the time the spins travel between sites, that is r has the value of the inverse 
of the band width. The resulting Green's function, which must be first smoothed out and 
normalised, interpolates between the VCA propagator at r = and the CPA propagator 
at t = oo. The density of states at finite r exhibits both the CPA sub-bands (or one CPA 
band, if U is small) and a VCA band. The bands are weighted according to the value of r. 
The CPA sub-bands are fixed at their positions at E — and E = U, while the VCA band, 
centred at E — elf, moves along with the concentration c. Equilibrium conditions with a 
polarised ground state are therefore found that are more stringent than the Stoner criterion 
for finite r. 

Artifacts of the r-CPA mixture include some non-regular features in the density of states 
for a range of parameters (intermediate U and concentration) , and peculiar behaviour of the 
m versus U curves at large r and small concentration. The r-CPA is however a successful 
method for the Hubbard model that is based on the better physical picture given by the 
combination of the CPA and VCA. It allows the treatment of the scattering by the CPA 
while also predicting ordered ground state for large enough U. In a next publication, it will 
be shown that the r-CPA can also be used to calculate two-particle Green's functions. In 
particular the dynamical susceptibility in the Hubbard model can be evaluated using this 
method. 
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